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Abstract 

We consider a finite dimensional damped second order system and 
obtain spectral inclusion theorems for the related quadratic eigenvalue 
problem. The inclusion sets are the 'quasi Cassini ovals' which may 
greatly outperform standard Gershgorin circles. As the unperturbed 
system we take a modally damped part of the system; this includes 
the known proportionally damped models, but may give much sharper 
estimates. These inclusions are then applied to derive some easily cal- 
culable sufficient conditions for the over damp edness of a given damped 
system. 

1 Introduction and preliminaries 

A damped linear system without gyroscopic forces is governed by the differ- 
ential equation 

Mx + Cx + Kx = f{t). (1) 

Here x = x{t) is an R"-valued function of time t G M; M, C, K are real 
symmetric matrices of order n. Typically M, K are positive definite whereas 
C is positive semidefinite. The physical meaning of these objects is 

x{t) position or displacement 
M mass 

C damping 

K stiffness 
f{t) external force 
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If in the homogeneous equation above we insert x(t) = e^^x, x constant, we 
obtain 

{\^M + \C + K)x = Q (2) 

which is called the quadratic eigenvalue problem, attached to ([1]), A is an 
eigenvalue and x a corresponding eigenvector. 

The quadratic eigenvalue problem may have poor spectral theory in spite 
of the hermiticity and positive (semi)definiteness of M.C, K. There always 
exists a non-singular matrix $ such that 

$^M$ = /, <I>^K$ = Vl = diag(u;^, . . . , cj^). (3) 

If the matrix $ can be chosen such that also 

D = <I>^C$ (4) 

is diagonal then the system is called modally damped. 

While (12]) is the standard spectral decomposition of a symmetric positive 
definite matrix pair, a simultaneous achieving of (jlj) is rather an exception 
being equivalent to the generalised commutativity property 

CK^^M = MK~^C. (5) 

However, as an approximation, modal damping is attractive since it is han- 
dled by the standard theory and numerics of Hermitian matrices. The aim 
of this paper is to assess modal approximations of general damped systems. 
More precisely, we will derive spectral inclusion theorems for eigenvalues 
where the unperturbed system is modally damped. There is some hierarchy 
among various modal approximations of a given damped system and we will 
investigate this issue as well. Our inclusion sets will not be circles, we will 
call them quasi Cassini ovals. We will show that our ovals outdo classical 
Gershgorin circles. A special case are overdamped systems the eigenvalues of 
which are particularly well behaved, there ovals reduce to intervals and inclu- 
sions of Wielandt-Hoffman type will be derived. Finally, we will derive new 
calculable sufficient conditions for the overdampedness of a given system. 

2 Modal approximation 

Some, rather rough, facts on the positioning of the eigenvalues are given 
in [3]. Further, more detailed, information is obtained by the perturbation 
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theory. A simplest thoroughly known system is the undamped one. Next to 
this lie the modally damped systems. 

A simplest eigenvalue inclusion for a general matrix A close to a matrix 
Ao is 

(r{A)CS, = {\: \\{A-Ao){Ao-\I)-'\\<l} (6) 
Obviously Si C ^2 with 

32 = {A: \\{A,-XI)-T'<\\{A-Ao)\\}. (7) 

This is valid for any matrices A, Aq. Using $, Q from we set 

so the quadratic eigenvalue equation ([2]) is equivalent to 

Ay = Xy. 



Here we have set 



A 








, ^0 = 





Q ' 


-n 


-D 


-n 






(8) 



(9) 



Hence 



A -A, 




-D 



The matrix Aq is skew-symmetric and therefore normal, so || {Aq — XI) 
dist(A, (t(Ao)) hence 



-iii-i 



32 = {A: dist(A,or(Ao) < UA-A, 



where 



II (^-^o) 



II^I 



1-^2 "^C'-Lj '^1 



x^Cx 



max ■ 



(10) 
(11) 



x'^Mx 

is the largest eigenvalue of the matrix pair C, M. We may say that here 'the 
size of the damping is measured relative to the mass'. 

Thus, the perturbed eigenvalues are contained in the union of the disks 
of radius \\D\\ around a{Ao). 
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Remark 2.1. In fact, cr{A) is also contained in the union of the disks 

{A : \XTii^j\ < Rj} (12) 

with 

n 

Rj = Y,\d,,\. (13) 

k=l 

(Replace the spectral norm in (jS]) by the norm || • ||i). 

The bounds obtained above are, in fact, too crude, since we have not 
taken into account the structure of the perturbation A — Aq which has a 
remarkable zero pattern. 

Instead of working with the matrix A we may turn back to the original 
quadratic eigenvalue problem in the representation in the form (see ([3]) and 

©) 

det{\'^I + \D + n^) =0. 



The inverse 



exists, if 



(A'/ + + XD{x^i + nYY^ 

\\D{X^I + n^)-'\\\X\ <1 (14) 
which is implied by 

\\ix'i + nr'\\\\D\\\x\ = . ,1 ''^'"hL • n < 1 (15) 

mmj(|A — tujj\\X + tujj\) 

Thus, 

a{A) C UjGiiujj, -iuoj, \\D\\), (16) 

where the set 

e(A+,A_,r) = {A : |A-A+||A- A_| < |A|r} (17) 

will be called quasi Cassini ovals with foci X± and extension r. This is in 
analogy with the standard Cassini ovals where on the right hand side instead 
of |A|r one has just r^. (The latter also appear in eigenvalue bounds in 
somewhat different context.) We note the obvious relation 

e(A+, A_, r) C e(A+, A_, r'), whenever r < r'. (18) 
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The quasi Cassini ovals are qualitatively similar to the standard ones; they 
can consist of one or two components; the latter case occurs when r is suf- 
ficiently small with respect to |A+ — A_|. In this case the ovals in f|T6l) are 
approximated by the disks 

\X±^u;,\<^^ (19) 

and this is one half of the bound in (fTOl) . ( fTTl) . 

Remark 2.2. (j[A) is also contained in the union of the ovals 

e{iujj,-iLJj, \\Q-^DQ-'^\\uj^). (20) 
Indeed, instead of inverting X^I + XD + il"^ invert X^il'^ + Afi"^Df2-^ + I. 
Remark 2.3. (J {A) is also contained in the union of the ovals 

G{iujj, —iujj, Rj) (21) 

and also 

G{iLjj, —iuj, pjUjJ) (22) 

with 

ft = E^- (23) 

The just considered undamped approximation was just a prelude to the 
main topic of this section, namely the modal approximation. The modally 
damped systems are so much simpler than the general ones that practitioners 
often substitute the true damping matrix by some kind of 'modal approxi- 
mation'. Most typical such approximations in use are of the form 

Cprop = aM + (3K (24) 

where are chosen in such be in some sense as close as 

possible to C, for instance, 

Tr [{C -aM- (3K)W{C - aM - ^K)] = min, (25) 

where W is some convenient positive definite weight matrix. This is a propor- 
tional approximation. In general such approximations may go quite astray 
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and yield thoroughly false predictions. We will now assess them in a more 
systematic way. 



A modal approximation to the system ([T]) is obtained by first representing 
it in modal coordinates by the matrices D, VL and then by replacing D by its 
diagonal part 

D° = diag(dn,...,d„„). (26) 

The off-diagonal part D' = D — is considered a perturbation. Again 
we can work in the phase space or with the original quadratic eigenvalue 
formulation. In the first case we can make perfect shuffling to obtain 



Ui 
—uji da 



A, 



A — (Aij), An - 

Aq = diag(Aii, . . .,Ann)- 





d. 



(27) 
(28) 



So, for n = 3 
























-UJi 


-dn 





-di2 





-di3 


A = 











UJ2 











-di2 


-UJ2 


-d22 





-d23 



















UJ3 







-di3 







-UJ3 


-d33 



Then 



ll(^-A/)-^| 



max II {A 

j 



A/)-| 



Even for 2 x 2-blocks any common norm of {Ajj — XI) ^ seems complicated 
to express in terms of disks or other simple regions, unless we diagonalise 



Xi 

xi_ 



xi 



(29) 



As is directly verified, 




^33 

2ui 
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with 

9, = 

Set 5* = diag(S'ii, . . . , Snn) and 

A' = S-^AS = A'o + A" 

then 

A'o = diag(Ai,...,A^), 
Ajf^ = Sj Ajf^Sfs, A = S AS 
Now the general perturbation bound (fTOi) . apphed to A'q,A", gives 

(t(A) C Uj-±{A : |A- Ail < (30) 

There is a related 'Gershgorin-type bound' 

a(A)CU,,±{A: \X - Xi\ < ^{S,)r,} (31) 

with 

n 

r,=J2\\djk\\. (32) 

To show this we replace the spectral norm || ■ || in (I6l) by the norm ||| • 
defined as 

:= max^^ ll^fcill 

k 

where the norms on the right hand side are spectral. Thus, ([6]) will hold, if 

max^||(A-Ao)fc,||||(%-A/)-i||<l 

Taking into account the equality 

ll(^-^o).,|| = { '^^^■[j It'j 

A G cr{A) implies 

r,>||(^,,-AJ)-'||> '°'°'l^-^||-l^-^-l' 



and this is ([21]) • 



Note that the bounds fl30l) and fl3Tl) are poor whenever the modal approx- 
imation is close to a critically damped eigenvalue. 

Better bounds are expected, if we work directly with the quadratic eigen- 
value equation. The inverse 

(A'/ + AZ}° + + \D\xH + AD° + ^^)-^)-^ 

exists, if 

\\D'{\^I + \D^ + n'')-^\\\\\<l (33) 

which is insured, if 

+ AD° + n^)-^\\\\D'\\\\\ = \\£3!^ _ < 1 (34) 

' ' min,(|A- A^+||A- Ail) ^ ' 

Thus, 

a(A)CU,e(Ai,Ai,||D'||). (35) 

These ovals will always have both foci either real or complex conjugate. If 
r = \\D'\\ is small with respect to [A-^ — Ai| = \J\<Pjj — 4:U!'j\ then either |A— A-'^l 

or |A — A-L| is small. In the first case the inequality |A — A-'^HA — A-L| < |A|r 
is approximated by 



A-Ai|<-f-^ = r ^ (36) 



and in the second 



djj > 2uj 



da < 2u!i 



Wnr I ^^^^ ' 
|A-Ai|< ' ' . =rl . (37) 

V JJ J 



This is again a union of disks. If djj ~ then their radius is ~ r/2. If 
djj ^ 2ujj i.e. A_ = A_|_ ^ —djj/2 the ovals look like a single circular disk. 
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For large djj the oval around the absolutely larger eigenvalue is ~ r (the 
same behaviour as with (1311) ) whereas the smaller eigenvalue has the diameter 
~ 2rLu'j/d'jj which is drastically better than fl3Tl) . 

In the same way as before the Gershgorin type estimate is obtained 

a{A) CU^eiX^, XI, Vj). (38) 

We have called D' a modal approximation to D because the matrix D 
is not uniquely determined by the input matrices M, C, K. Different choices 
of the transformation matrix $ give rise to different modal approximations 
D' but the differences between them are mostly non-essential. To be more 
precise, let $ and $ both satisfy ([3]). Then 

K = = ^-'^n'^^-^ 

implies that f/ = is an orthogonal matrix which commutes with 

Q = dmg{LJiIn^, . . . ,UjJns), UJi<---<Us. (39) 

Hence 

t/° = diag(?7n,...,f/.s), 
where each Ujj is an orthogonal matrix of order rij from Now, 

D = ^'^C^ = U^^'^C^U = U^DU, (40) 

A, = U^D,,U, (41) 

and hence 

D' = U^D'U. (42) 

Now, if the undamped frequencies are all simple, then U is diagonal and the 
estimates (fTTj) or (fT6l) - (fT7|) remain unaffected by this change of coordinates. 
Otherwise we replace diag((iii, . . . , dnn) by 

A = diag(Dn,...,As) (43) 

where commutes with Q. In fact, a general definition of a modal approx- 
imation is that it 

1. is block-diagonal and 
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2. commutes with Q. 



The modal approximation with the coarsest possible partition — this is the 
one whose block dimensions equal the multiplicities in Q — is called a max- 
imal modal approximation. Accordingly, we say that C° = ^~^D^^^^ is a 
modal approximation to C (and also M, C^, K to M, C, K). 

Proposition 2.4. Each modal approximation to C is of the form 

s 

C^ = Y^P*CPk (44) 

k=l 

where Pi, . . . .Ps is an M -orthogonal decomposition of the identity (that is 
P^ = MPkM~^) and Pk commute with the matrix 



VM-^K = M-i/VM-i/2i^M-i/2M^/2 
Proof. Use the formula 



D' = J2P'kDP', 



k=l 

with 

P° = diag(0 . . . , 4,, . . . , 0, D = $^C$, D° = $^C°<I> 

and set Pk = $Pfc°$-^ Q.E.D. 

It is obvious that the maximal approximation is the best among all modal 
approximations in the sense that 

\\D - D^We < \\D - D'^We, (45) 

where 

^° = diag(Dn,...,4.) (46) 

and D = (Dij) is any block partition of D which is finer than that in (H3|) . 
We will now prove that the inequality fl4S]) is valid for the spectral norm also. 
We shall need the following 

Proposition 2.5. Let H = (Hij) be any partitioned Hermitian matrix such 
that the diagonal blocks Ha are square. Set 

if° = diag(ifn, . . . , i^ss), H' = H- . 
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Then 

Xk{H) - UH) < \,{H') < \,{H) - \^iH) (47) 

where Afc(-) denotes the non- decreasing sequence of the eigenvalues of any 
Hermitian matrix. 

Proof. By the monotonicity property (Wielandt's theorem) we have 

\k{H) — maxmax(T(-f/jj) < Xk{H') < Xk{H) — min min a {Hjj). 

j j 

By the interlacing property, 

Xi{H) < a{Hj,) < Xn{H). 
Together we obtain (gT]). Q.E.D. 

From p71) some simpler estimates immediately follow: 

ll^^'ll < Xn{H) - Xi{H) =: spread(if) (48) 
and, if H is positive (or negative) semidefinite 

\\H'\\ < \\H\\. (49) 



Now fl45l) for the spectral norm immediately follows from fl49|) . So, a best 
bound in fl35l) is obtained, if = D — D' is a maximal modal approximation. 

Proposition 2.6. Any modal approximation is better than any proportional 
one. 

Proof. With 

Dprop = C(I + (3Q 

we have 

\iD - Dprophl > \iD - D%\ = \D: 

which implies 



\D — DpfopW ^ ll-D I 



Q.E.D. 



If is block diagonal and the corresponding D' = D — is inserted in 
fl5S]) then the values djj from fl2^ should be replaced by the corresponding 
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eigenvalues of the diagonal blocks Djj. But in this case we can further 
transform Q and D by a unitary similarity 

U = dmg{U,,...,Us) 

such that each of the blocks Djj becomes diagonal {Q stays unchanged). 
With this stipulation we may retain the formula fl35l) unaltered. This shows 
that taking just the diagonal part D° of D covers, in fact, all possible modal 
approximations, when $ varies over all matrices performing ([3]). 

Similar extension can be made with the bound ( l38l) but then no improve- 
ments in general can be guaranteed although they are more likely than not. 

By the usual continuity argument it is seen that the number of the eigen- 
values in each component of UiC{X\, X^_,ri) is twice the number of involved 
diagonals. In particular, if we have the maximal number of 2n components, 
then each of them contains exactly one eigenvalue. 

A strengthening in the sense of Brauer is possible as well. We will show 
that the spectrum is contained in the union of double ovals, defined as 

D(A^+,A^_,A^,Al,r,r,) = 

{A : |A-A^||A- A^IIA- A^IIA- All < rpr,|Ap}, (50) 

where the union is taken over all pairs p q and A^ are the solutions of 
A^ + dppX + ujp = and similarly for A?^. The proof just mimics the standard 
Brauer's one. The quadratic eigenvalue problem is written as 

n 

(A^ + Xdpp + uj^)xi = -A ^ dijXj, (51) 

i=i 
n 

(A^ + Xdii + ujf)xi = —A dijXj, (52) 

3 = 1 

where \xp\ > \ Xq I dbTQ the two absolutely largest components of x. If Xq = 
then Xj = for all j ^ p and trivially A G 2)(A!J_, A^, A^, Ai, rpr^). If Xg 7^ 
then multiplying the equalities (1^ and (13^ yields 

|A - A^IIA - A^ ||A - A^IIA - X'L\\xp\\xg\ < 



12 



3=1 k=l 



Because in the double sum above there is no term with j = k = p we always 
have < |xp||xq|, hence the said sum is bounded by 



I I I "^p 1 1 "^g I ^ ^ I ^pj I ^ ^ I dqk I • 



j=i fc=i 



Thus, our inclusion is proved. As it is immediately seen, the union of all 
double ovals is contained in the union of all quasi Cassini ovals. 



The simplicity of the modal approximation suggests to try to extend it 
to as many systems as possible. A close candidate for such extension is any 
system with tightly clustered undamped frequencies, that is, Q is close to an 

from (!39l) . Starting again with 

(A^j + + + {XD' + z) + {x^i + + 

with Z = ~ (Q^)"^ we immediately obtain 

a{A)CU,e{Xi,Xl,\\D%\\Z\\). (53) 

where the set 

e(A+,A_,r,g) = {A : |A- A+||A- A_| < |A|r + g} (54) 
will be called modified Cassini ovals with foci A± and extensions r, q. 

Remark 2.7. The basis of any modal approximation is the diagonalisation 
of the matrix pair M, K. An analogous procedure with similar results can 
be performed by diagonalising the pair M, K or C, K. 



3 Modal approximation and overdampedness 

If the systems in the previous section are all overdamped then estimates are 
greatly simplified as ovals become just intervals. But before going into this a 
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■12 « 411 111 4H -1.1! -1.1 -11 U W it -1 4! 41 4J 4.S 45 44 43 42 41 

Figure 1: Ovals for u; = 1; d = 0.1, 1; r = 0.3 
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-If -W -12 -I M M -M 1 ■ ■! -li -I * 1 ■ -! -li -1 

Figure 2: Ovals for u; = 1; d = 1.7, 2.3, 2.2; r = 0.3, 0.3, 0.1 
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more elementary — and more important — question arises: Can the modal 
approximation help to decide the overdampedness of a given system? 

We begin with some obvious facts the proofs of which are left to the 
reader. 

Proposition 3.1. If the system M, C, K is overdamped, then the same is 
true of the projected system 

M' = X*MX, C' = X*CX, K' = X*KX (55) 

where X is any injective matrix. Moreover, the definiteness interval of the 
former is contained in the one of the latter. 

Proposition 3.2. Let 



M = 


diag(Mii, . 


• . , Mss) 


C = 


■■ diag(Cii, . 


■ ■ ) Css) 




diag(Kii, . 


...Kss). 



Then the system M, C, K is overdamped, if and only if each of Mjj, Cjj, Kjj 
is overdamped and their definiteness intervals have a non trivial intersection 
(which is then the definiteness interval of M, C,K) 

Corollary 3.3. // the system M, C, K is overdamped, then the same is true 
of any of its modal approximations. 

Obviously, if a maximal modal approximation is overdamped, then so are 
all others. 

In the following we shall need some well known sufficient conditions for 
negative definiteness of a general Hermitian matrix A — {aij)] these are: 

%i < (56) 

for all j and either 

\\A - diag(aii, . . .,ann)\\ < -maxa^.,- (57) 

3 

(norm-diagonal dominance) or 

n 

\akj\ < —djj for all j (58) 

fc=i 

(Gershgorin-diagonal dominance). 
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Theorem 3.4. Let Q, D, rj be from ^3^, respectively and 

D° = diag(dn ,...,dnn), D' = D-D^. 

Let either 

Aj = {djj - II 11)2 - Alu] > for all j (59) 

and 

-djj + IID'II - a/A~ . -djj + \\D'\\ + ^/A'j 



II II \ J ^ • JJ ' II II ' V J /nn\ 

p_ := max < mm =: (oUj 

j 2 j 2 

or 

Aj = {djj -Tjf - Auj] > for all j (61) 

and 

-djj + Tj-Jkj ^ -djj + rj + JAj 

p_ := max < mm =: (62) 

j 2 j 2 

Then the system M,C,K is overdamped. Moreover, the interval (p_,p+), 
{p-,p+), respectively, is contained in the definiteness interval of M,C, K . 

Proof. Let p^ < fi < p-^-. The negative definiteness of 

will be insured by norm-diagonal dominance, if 

— /i||D'|| < — — fidjj — uj'j for aU j, 
that is, if /i hes between the roots of the quadratic equation 

/i^ + ^i{djj - ||D'||) + = for all j 



and this is insured by (1591) and ( l60l) . The conditions ( pTj) and ( 1621) are treated 
analogously. Q.E.D. 

We are now prepared to adapt the spectral inclusion bounds from the 
previous section to overdamped systems. Recall that in this case the defi- 
niteness interval divides the 2n eigenvalues into two groups: J-negative and 
J-positive. 
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Theorem 3.5. // ^5M) o,nd / fgOj) hold then the J -negative/ J -positive eigen- 
values are contained in 

U, U,(/ii_, (63) 

respectively, with 



-d,,-\\D'\\±J{d,, + \\D'\\f-Auj^^ 
f^i. = (64) 



-dj, + \\D'\\±Jid,,-\\D'\\r-Au;] 
/^+- = • (65) 

An analogous statement holds, if l[61\) and [6^) hold and , is replaced 
by fji!'++ , where in ([g^|fg3]) ||-D'|| is replaced by rj. 



Proof. All spectra are real, so we have to find the intersection of C(A'?,_, X'L,r) 
with the real line the foci A'^, A'L from (!29|) being also real. This intersection 
will be a union of two intervals. For A < A''_ and also for A > A'^ the j-th 
ovals are given by 

(Ai-A)(A^+-A) < -Ar 

i.e. 

A' - (Xi + Ai - r)A + A^Ai < 

where X'_^_ + Xi_ = —djj and X-'^X'L = cuj. Thus, the left and the right boundary 
point of the real ovals are ++ . 

For A-L < A < A-^ the ovals will not contain A, if 

(A- Ai)(A^+-A) < -Ar 

i.e. 

A^ + {djj - r)A + cj| < 

with the solution 

< A < 

Now take r = ||-D'||. The same argument goes with r = rj. Q.E.D. 
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Note the inequality 

(66) 

for all j, k. 

Monotonicity-based bounds. As it is known for symmetric matrices 
monotonicity-based bounds for the eigenvalues (Wielandt-Hoffmann bounds 
for a single matrix) have an important advantage over Gershgorin-type bounds: 
While the latter are merely inclusions, that is, the eigenvalue is contained 
in a union of intervals the former tell more: there each interval contains 'its 
own eigenvalue', even if it intersects other intervals. 

In this section we will derive bounds of this kind for overdamped systems. 
A basic fact is the following theorem 

Theorem 3.6. With overdamped systems the eigenvalues go asunder under 
growing viscosity. More precisely, Let 

A-_„<---<Ar<A+<---<A+ <0 

be the eigenvalues of an overdamped system M, C,K. If M, C, K is more 
viscous that is, M < M, C > C,K < K in the sens of forms then its 
corresponding eigenvalues A^ satisfy 

k<K, A+<A+ (67) 

A possible way to prove this theorem is to use the Dufiin's minimax 
principle pj, moreover, the following formulae hold 

A^ = minmaxp+(a;), A^ = maxminp_(x). (68) 

Sk xeSk Sk x&Sk 

where is any fc-dimensional subspace. Now the proof of Theorem 13.61 is 
immediate, if we observe that 

p±{x)^p±{x) (69) 

for any x. 

As a natural relative bound for the system matrices we assume 
\x'^6Mx\ < ex^Mx, \x^6Cx\ < ex^Cx, \x^6Kx\ < ex^Kx, (70) 
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with 

SM^M -M, dC^C -C, SH^K -K, e < 1. (71) 

We suppose that the system M, C, K is overdamped and modally damped. 
One readily sees that the overdampedness of the perturbed system M , (7, K 
is insured, if 

d—1 , . x^Cx ,„^, 

e<- , ci = mm — =. (72) 

So, the following three overdamped systems 

(1 + e)M, (1 - e)C, (1 + e)i^; M,C,k- (1 - e)M, (1 + e)C, (1 - e)i^ 

are ordered in growing viscosity. The first and the last system are over- 
damped and also modally damped and their eigenvalues are known and given 

by 

l + e 



respectively, where 



K iv) = , 

are the eigenvalues of the system M, rjC, K. We suppose that the unper- 
turbed eigenvalues = A^(l) are ordered as 

A" < • • • < A7 < At < • • • < A+. 

n — — 1 1 — — n 

By the monotonicity property the corresponding eigenvalues are bounded as 
\ti\^\<'>t<\t{\^\, (73) 



-l + e, 

where X^{r]) are obtained by permuting X'^{r)) such that 

A^(^)<---<A+(r;) 

for all T) > 0. It is clear that each A^(77) is still non-decreasing in 77. An 
analogous bound holds for as well. 
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